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Abstract. A new class of semi-analytically solvable MHD a 2 — dynamos is found 
based on a global diagonalization of the matrix part of the dynamo differential 
operator. Close parallels to SUSY QM are used to relate these models to the Dirac 
equation and to extract non-numerical information about the dynamo spectrum. 
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The magnetic fields of planets, stars and galaxies are maintained by dynamo 
effects in electrically conducting fluids or plasmas. A crucial role in the qualitative 
understanding of the basic dynamics of dynamos play various toy models which allow 
for a semi-analytical study. In this respect, the spherically symmetric a 2 — dynamo 
[HE] plays a role similar to the harmonic oscillator in Quantum Mechanics (QM). But 
even this model can be analytically described only in very few cases — like, e.g., for 
constant a— profiles (The a— profile a(r) acts as an effective potential in the dynamo 
operator.) or in the case of idealized boundary conditions [3] (mimicking a perfect 
external conductor [4j or the limit of high angular mode numbers I 3> 1 [5]). 

Here, we are going to present another class of highly simplifying a 2 — dynamo 
models based on a global diagonalizability of the matrix part of the dynamo differential 
operator. For this purpose we relax the very rigid boundary conditions (BCs) at 
the surface of the dynamo maintaining fluid and replace it by a combination of a 
strongly localized a embedded in a conducting surrounding [5] and Dirichlet BCs at 
infinity. This allows us to use not only the Krein space symmetry properties of the 
operator (closely related to VT— symmetric QM [9]) but also to uncover deep 

relations to super-symmetric (SUSY) QM (for a recent review see [10]) and to map the 
dynamo eigenvalue problem into a set of two 2x2 stationary Dirac equations. Similar 
to [B], we will find a close relation between overcritical dynamo states and bound 
states in an associated QM model. For a certain parameter value the diagonalization 
technique breaks down and we develop a perturbative approach to describe the system 
behavior in the vicinity of the corresponding operator Jordan structure. The spectral 
reality properties of overcritical dynamo states are discussed. Finally, we interpret our 
findings in terms of a special link between the radial components of the electromotive 
force and the induced currents. 
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Subject of our analysis is the eigenvalue problem of a spherically symmetric 
a 2 — dynamo 2 in its simplified and unitarily re-scaled form [71(5] 

(d r Kd r + M - XI) $(r) = , $(0) = , $(oo) = 0. (1) 

The matrix structure is encoded in 

K =I-aa-, (2) 

M = - K l -^±^ + aa+ , (3) 

where / is the 2x2 unit matrix and a± denote the nilpotent matrices o± — ^{o- x ±io~ y ) 
with a x ,y,z being the usual Pauli matrices. The system (fTJ) describes the coupled 
poloidal and toroidal components of the /—modes of the magnetic field in a mean- 
field a 2 — dynamo model with helical turbulence function (a— profile) a(r). We 
assume this a— profile real-valued, bounded, sufficiently smooth and exponentially 
decreasing for large r> 1. The latter assumption allows us to consider the a— profile 
approximatively as localized and at the same time to relax the otherwise rigidly 
imposed physical BCs at a given fluid/plasma surfaccFil (set, e.g., at r = 1) replacing 
them by Dirichlet BCs at r — > oo. Such an approach will make underlying structural 
links to quantum mechanical setups transparent and will provide the Krein-space self- 
adjointness of the eigenvalue problem |T]) similar to that of models with idealized 
(Dirichlet BCs) at fixed r = 1 (see Refs. [3 Hi]). 

Our first goal is in diagonalizing the matrix structure of the 2x2 matrix differential 
operator in {!]). For this purpose we use a two-step procedure, which consists in a 
replacement of the dependent variable by a Kummer-Liouville type transformation 

$(r) = P _1 (r)E(r) (4) 

to remove acr_ from K in the derivative term and an afterwards performed coordinate- 
independent (global) similarity transformation to diagonalize the remaining matrix 
potential. The matrix P can be found to be a square root of K, 

,2 „ / 1 



where 



K = P\ P=i* I ), (5) 
and yields the following equation for S(r): 

s » (r) _ja+i) s(r)+7(r)s(r)=0i (6) 

/ ia 2 (r)-A a(r) \ 

V(r) = 2 • (7) 

\ ±a"{r) + \a 3 {r)-a{r)\ \a 2 {r) - X J 

The equation system (j^J can be globally decoupled provided the eigenvectors of the 
matrix (O are r-independent. It is not difficult to see that this is possible if the 
function a(r) satisfies the equation 

a"(r) + ^a 3 (r)-a 2 a(r)=0 (8) 

with a an arbitrary real constant. The general solution to this equation is expressed in 
terms of an elliptic integral which under the additional requirement \a(r — > oo)| — > 
reduces to 

/ \ 2a 

a r - —rr- ( V] (9) 

cosh|a(r — r^)\ 



1 ) The physics of MHD dynamos is discussed, e.g., in [TJ [2] . 
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with ar Q as an integration constant. For a— profiles © the matrix part of the dynamo 
operator can be diagonalized — except for the special case A = \a? when it is similar 
to a 2 x 2 Jordan block. Under the diagonalization V i— » U~ 1 VU Eq. © splits into the 
following decoupled pair of differential equations (DEs) for the components (F + , F_) 
of the vector If- 1 -. = (F + ,F-) T : 



-d 2 + - -a 2 T ea 

r l 2 

,V2„„„ / 1 9 x\V 2 



F± = -\F± (10) 



with e = (\a 2 — A) " or e = — {^a 2 — A) and the diagonalizing matrix U given by 

"-{]!.)■ ("J 

We notice that in equations © and (|10[) the parameter a is inessential and can be 
eliminated by re-scaling e — as, r = x/a, A = a 2 X, a = aa. Further on, we will 
work in "a units" what is equivalent to setting a — 1 and identifying r — x, e = i, 
A = A, a = a. Apart from PU|) we will also use reshaped versions of these equations 
(obtained by substitution of A = i — e 2 ) which take the form of quadratic pencils in 
the auxiliary spectral parameter e, 

l-d 2 x + l -V±V- 1 -a 2 + 1 - T ea-e 2 }F ± = 1 a = — -| (12) 

x £ 2 2 cosh(x — xo) 

supplemented by the Dirichlet BCs F±(x = 0) = F±(x = oo) = 0. 

In the special case e = 0, i.e. for A = Xj := |, the diagonalization matrix J7 (see 
(llip ) becomes singular and the system © assumes the upper triangular (Jordan- type) 
form 

d 2 -v -v, \/sr 

9^ — / I S 



(13) 



with 3q, Si as components of the vector S = (Si, So) and the potential terms given 
by Vo = l(l + l)x~ 2 -Ua 2 -1), V 1 



-a. 



We start our investigation of Eqs. (flTJ)) . JT2T) by noticing that for I ^ they are 
DEs of non-Fuchsian type and therefore their solutions cannot be expressed in terms 
of ordinary special functions. A certain simplification occurs for the monopole case 
1 = 0. Then (fT0|) . (fT2|) constitute a particular type of Heun's equations having three 
different finite regular singularities (see e.g. [11]) and solutions expressible in terms of 
Heun's functions. A corresponding analysis will be presented elsewhere. 

Here we are going to use the fact that Eqs. (P~0|) . (I12p are closely related to the 
exactly solvable stationary Schrodinger equation 

H 1( f> = E4>, Hl = -d 2 x -^a 2 (14) 

with Hi well known as superpartnejf}] of the trivial Hamiltonian Hq — —d 2 

LHq = H\L, L = —d x + w , w = d x hiu, u — cosh(a; — xo),(15) 
L f L = H + 1, LL^ =H! + l, (H - E f )u = 0, E f = -1 (16) 

so that for E > Ef = — 1 the solutions on the halfline x £ [0, oo) are simply given as 
(/) ± =Le ±KX = [Tn + t&nh{x-xo)}e ±KX , k = V^E . (17) 



2 ) The potential — « 2 /2 = —2/ cosh 2 (a; — xo) is a xo— shifted modified Poschl- Teller [12] and reduced 
Rosen-Morse potential 1131 and also known as "one-soliton potential well" in KdV theory. 
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Imposing the Dirichlet BC 4>(x — 0) = and exponential decay on these solutions 
makes Hi an essentially selfadjoint operator with a single discrete level E(xo) = 
— tanh 2 (a;o) € (—1,0) for x > and with the halfline E > as its continuous 
spectrum. The additional dynamo-related constraint <p(x — > oo) — > selects then the 
bound state (BS) as relevant solution. With this information at hand on the spectrum 
of the exactly solvable Schrodinger equation (jT4j) , we are now well prepared to present 
a qualitative discussion of the interrelated spectra of the dynamo eigenvalue problems 

CD). 

A first piece of information can be extracted from (|10p by neglecting for a moment 
the potential terms 1(1 + l)/x 2 and +ea . In this case (fTOf structurally coincides with 
(|14p and we can identify E = —A. This means that due to the physical constraint 
<j)(x — > oo) — * and its implication i? 6 (—1,0) the model necessarily describes 
overcritical dynamo regimes 3 )! A > 0. 

In order to extract further information, we proceed with the familiar QM model 
and extend Hi by the centrifugal potential 1(1+1) /x 2 . This potential acts as a repulsive 
barrier in the vicinity of the origin x = and for small Xo it overcompensates the effect 
of the attractive one-soliton potential well —a 2 /2 = —2/ cosh 2 (x — xq) with center at 
x = xq. As result, no BS can exist for small xq. When xo is increased beyond a certain 
/—dependent critical value the effect of the repulsive barrier will become sufficiently 
weak and the BS level will re-appear from the lower boundary E — of the continuous 
spectrum and move down toward the lower boundary E — > Ef = — 1 of the BS band. 
The situation is illustrated in Fig. [Tk showing the sign-inverted picture for A = — E. 




Figure 1. Spectra A(xo) for the reduced system [—d 2 + 1(1 + l)/x 2 — 
a 2 /2]F± = — \F± (a) and for the complete problem (|10[1 (b) in case of 
angular mode numbers I — 0, 1,2,3. For numerical reasons the Dirichlet 
BC has been imposed at the large distance x = 100. 

It remains to clarify the role of the energy dependent part +sa = +(\ — A) 1 ' 2 a 
of the potential. Due to the specific square-root coupling this term might, in general, 
become complex-valued. Hence, the main question to answer is whether the spectrum 
will remain purely real (as the numerical results in Fig. [TJd indicate) or whether and 
under which conditions it might become complex. 

A partial clarification can be achieved by transforming the pencil equations ()12|) 
into equivalent Dirac equations. For this purpose we use a more general and slightly 
reshaped version of the SUSY factorization technique (fT6|) (cf. p3]), factorizing instead 

3 ) A dynamo in its kinematic regime is called overcritical when it has a positive growth rate 
Re(A) > 0. This is in contrast to so called undercritical regimes Re(A) < with decaying (dissipating) 
magnetic field. (See e.g. [HIS].) 
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of Hi the Hamiltonian 

Hv = -%+ l -£±V-±c? + \=ltL, (18) 

L = —d x + w, L'=d x +w, w — u/u, (19) 

which has a continuous spectrum for E > 1/2 and a BS at some E < 1/2. The 



factorization (|T8j) allows us to rewrite the pencils (|T2|1 as 

[L + LT£a-e 2 ]F± = 0. (20) 
Denoting now _F± =: -0i and assuming for e =/= 

L^i =: (21) 

equation ([2U|l can be expressed as 

(±a + e)Vi =0. (22) 

From the explicit form (flT))) of the operators L and L + we see that equations (|2ip and 
(|22p are nothing but one-dimensional Dirac systems in their general representation 
(see e.g. [15]) 

= e* , i? = jd x + V (23) 



wit^ 



-($;)■ ^-(-i ;)■"-(- (24) 

For the factorization (fl~8|) to work, the function should be a solution of the 

equation H2,iu(x) — 0, i.e. an eigenfunction at the factorization energy Ef — 0. At 
the same time it should be nodeless on the positive semiaxis providing in this way a 
superpotential w(x) which is regular. According to an implication of the "oscillation" 
theorem (see, e.g., jT^]) a real eigenfunction of the operator (fTSj) can be nodeless only 
if its eigenvalue is located below the ground state level if the latter exists or below 
the lower bound of the continuous spectrum otherwise. Applied to our configuration 
with Ef = this means that u is nodeless as long as the BS (ground state) energy 
E(xo) > Ef. Comparison of (fl~8"|) with (fT2")) shows that Ef = coincides with the 
energy of the Jordan configuration (fT3")) at e — so that for xq < xj the solution u is 
nodeless and therefore w continuous. Together with the Dirichlet BC ipi(x = 0) = 
and the easily verified boundedness of \tp2(% = 0)| = |e _1 iV'iU=o < C < oo for 
£ / this makes an adapted version of Theorem III. 7.1 from 15j applicable: If the 
coefficients of a Dirac system are continuous functions in any finite interval of the 
positive semi-axis then the Dirac operator H defined by the differential system (f2"3"]l 
and the boundary condition ipi (0) cos 5 + "02 (0) sin i5 = (with S an arbitrary real 
number) is selfadjoint provided its domain of definition is ipi,ip2 G £ 2 (0, oo) with 
Lipi £ £ 2 (0,oo) , L+-02 — ottpi £ £ 2 (0, oo) . This means that the Dirac system (|2U| 
supplemented with the Dirichlet BC at the origin, and with it the dynamo problem, 
has a purely real spectrum for xq < xj. 

In case of xq > xj any real-valued solution to equation H2 } iu — has a node on 
the positive semi- axis so that the superpotential w(x) has a pole at some x > and 



4 ) With the help of a gauge transformation ^ = AO (see e.g. I15p which does not affect the zero 
boundary conditions, the Dirac equations J23H can be easily transformed into their canonical forms 
corresponding to a combination of scalar and pseudoscalar fields and describing the movement of 
massless particles. 
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the theorem is no longer applicable. A circumvention of this problem might consist 
in a complex linear combination u(x) = u^(x) + icuP"'{x), c € K. of two real-valued 
linearly independent solutions uW and uW of equation H 2 ju — 0. Such a function u 
has no nodes on the positive semi-axis, but the superpotential w becomes a complex- 
valued function. The spectrum of such a Dirac operator needs a special analysis which 
will be presented elsewhere. 

Here we proceed with a general qualitative analysis. We start from the observation 
that for small real A < 1/2 the square root e is real. Introducing an auxiliary parameter 
b £ K and replacing the potential term ^fea in (TTU)) by ba the quadratic pencils ([12")) 
in e reduce to an auxiliary linear eigenvalue problem in a b— dependent X(xo,b) with 

1/2 

additional constraints 6 = =F£ = =f(1/2 — A) . The existence problem of the pencil 
solutions can be easily studied graphically in the (b, A)— plane. A solution exists if the 
plot of the numerically obtained \(xa,b) for given xo has an intersection point with 
the graphics of one of the constraints b = 6(A). The corresponding analysis shows 
that for the branch e > one has a BS with (F + ^ 0, F_ = 0) for xq < x.j and 
with (F + = ^ 0) for xq > xj. Due to the invariance of (TIT))) and (fT2"|) under 

the sign (branch) change (s,F±) t— > (— e, F T ) this can be re-interpreted as a single 
solution (F + ^ 0, F_ = 0) with e(x a < xj) > and e(x > xj) < 0. The latter 
interpretation is confirmed by a direct numerical analysis (see Fig. [TJd) of Eqs. (fTQ|) 
(cross-checked by numerics on the original matrix-operator problem |T])) and by the 
graphics of e(xo) in Fig. [2l Although a rigorous and complete analytical proof of the 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 



Figure 2. Spectrum in terms of e(a;o) for angular mode numbers I — 
0,1,2,3. 

result for all io £ I, I > is still missing, we now present a perturbative treatment 
of an / = model in the vicinity of the non-diagonalizable operator-configuration (|13[) 
at e = which explains the numerical results at least locally. 
For I — the Jordan chain equations (fT5)) have the form 

[-dl - a 2 /2 + 1/2] So = 0, [-Sg - a 2 /2 + 1/2] Si = aS (25) 

and due to their relation to the eigenvalue problem (fT4| . (Hi — E)cf> = 0, E = —1/2 
they are exactly solvable. From (JT7J) we immediately read off that for BSs 

S = C o 0- , K = 2~ 1/2 , C = const . (26) 

At the same time we find the value of the parameter xq — xj for which the Jordan 
structure occurs: E(xj) = —1/2 = — tanh (xj) yields xj = arctanh (2 -1 / 2 ). The 
second (associated) BS component Si is easily obtained by standard techniques for 




2 4 6 8 10 
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inhomogeneous ODEs (see, e.g., [T7] ) 



dx' — d>- I ach+cb-dx' 



(27) 



with W denoting the constant Wronskian W = W{4>+,4>-) = — 2~ 1 / 2 . In the limit 
x — ► oo the product </> + j Q a(j) 2 _dx' diverges so that the BS condition Si(oo) = implies 
Co = and the full Jordan chain solution reads simply S = (Si,So) T = (Ci<fi-,0) T . 
Comparing the chain equations (|25[) at So = with the pencil equations (|12[) . i.e. 
with 

- i« 2 + i T ea - e 2 ]F± =0, (28) 



shows that for BSs Eqs. (|25[) can be interpreted as effective smooth limit of (|28|) 
at e — > 0. Hence, a BS perturbation theory can be constructed simply on the 
decoupled scalar Eqs. alone. As small perturbation parameter we choose the 

Xq— distance S — x — xj from the Jordan configuration and expand a — aj—a'j5+. . ., 
e = ei<5 + e2<5 2 + . . ., i<± = Si + x±<5 + . . ., where aj = 2/ cosh(x — xj) and ei, 62, 
. . . are coefficients to be defined from the perturbation scheme. This yields the defining 
equation for the first-order corrections x± 

[-a 2 - a 2 /2 + 1/2] x± = -ffiSi, gx ~aj[a'jTei] (29) 

with solutions 

Ci " 



X± - C ± 0_ , w 



(f>- I g!<p + (f)-dx' - 4>+ gi4>-dx' 



(30) 



The BS condition x±( x ~ * 00 ) ~~ * can only be fulfilled if J °° gi4> 2 _dx' = 0. Hence, it 
fixes the parameter 

J aj<j)±dx' 2 

and with it e = ^5/2. Knowing that for fixed Xq the spectral parameter e is the 
same in the equations for F + and F_ we conclude that e — =F<5 /2 can be valid only 
for one of the signs and that therefore it acts as a selection rule. Full compatibility 
with the numerical results and Fig. [2] is established by choosing the F+— related BS 
for e > and xq < xj, i.e. for 5 < 0, so that e = —5/2 holds for all sufficiently 
small 5 G (—01,02) and provides a smooth connection between the e > and e < 
branches. At the same time it excludes a BS for the solution F_ . 

Inspection of the recurrence algorithm for the higher order corrections shows 
that at each order the highest-order coefficient enters its defining equation only 
linearly so that no square roots are involved which could produce complex-valued 
contributions. Together with the reality of all other ingredients (<^±, ctj, a'j, . . .) of 
these recurrence equations we conclude that Im(e/-) = 0, Vfc € Z + and no complex- 
valued BS-e can emerge from an BS within the convergence region of the series 
e = J2T=i e kS k - In this way we found an argumentation complementary to the Dirac 
equation based technique for xq < xj. Another argument explaining the reality of the 
BS eigenvalue follows from the fact that the dynamo operator (fT]) with Dirichlct BCs is 
necessarily self-adjoint in a Krein space [3j [Tl [8] , and, hence, a spectral real-to-complex 
transition requires two spectral branches of different Krein-space type to coalesce at 
some point in parameter space. Once locally only a single BS exists there is also no 
chance for a BS-related spectral phase transition to complex eigenvalues within the 
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convergence region of e = 53fcLi e k$ k ■ The question of whether complex-conjugate BS 
eigenvalue pairs might split off from the continuum remains still open. 

Finally we interpret the obtained solution behavior F + ^ 0, F_ = in terms 
of a special link between the magnetic field components of the dynamo. According 
to [1] [2] (cf. also the appendix of Ref. [7]) the poloidal and toroidal components 
of the Ith angular and nth radial modes of the multipole expanded fields Bp and 
B^ n) are given by B^' n) = -V x (r x V) if' n) (r,0), B { t l ' n) = -r x Vif' n) (r, 6) 
with scalar functions F^\r,9) — r _1 $^' 2 ra ^ (r)F i °(0) built from spherical harmonics 
Y?(6) = y/(2l + 1)/(4tt)P ; (cos6>) and the solution components $^ 2 ™ ) of problem (!]). 
In our case of only one BS solution the n-dependence reduces to a single term <&i 2 . 
Expressing the components F± with the help of Eqs. @ and (jTTJ) in terms of $i 



^-f^T^^l (32) 



r 

F±=± 2iL- V2 
and using the inverted relations [T] between the original magnetic field strength 
anc j ^ ne corresponding scalar functions f[ 1 ^\ r • B^'") = — 1(1 + l)Fi' n \ 

r ■ (V x B( ( '")) = -1(1 + l)F^ n) , we arrive for F_ = after multiplying §^ with 
-1(1 + 1)Y,° at 

- 1(1 + m F- = ~y £ r ■ [AM - (| + e ) B^ 1 )] ^ . (33) 

(The relation B = /iH and one of the Maxwell equations, V x H = j, have been used.) 
Eq. (|33[) has to be interpreted as special link between the induced current j and the 
spectrally shifted electromotive force aB in the present dynamo model. However, it 
can be shown that the field is not a Beltrami field VxB / /3B. 
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